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ABSTRACT 

The global structure of optically thin hot accretion disks with radial advection included has 
been investigated. We solve the full energy conservation equation explicitly and construct 
the radial structure of the disk. It is found that advection is a real cooling process and that 
there are two solutions co-exist for a given mass accretion rate less than a critical limit. One 
is fully advection cooling dominated and the other is dominated by local radiative cooling. 
The advection dominated accretion disks are hotter than the local cooling dominated disks; 
they are most probably in the two-temperature regime and effects such as electron-positron 
pair production and annihilation may need to be considered to study the microphysics of 
the hot plasma. However, the global disk structure will not be much affected by the local 
radiative process. 
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1. INTRODUCTION 

In standard accretion disk models it is assumed that the disk is axisymmetric, non 
self-gravitating, geometrically thin and rotating at a strictly Keplerian law. The viscosity 
which is responsible for both the energy generation and the angular momentum transfer 
is described with a phenomenological a-model prescription (Shakura & Sunyaev 1973). 
Futhermore, the accretion flow is cooled by purely local radiative cooling process which 
balances the viscous heating exactly. The last assumption is justified if the global heat 
transport (radial advection) is relatively negligible. However, in optically thick slim ac- 
cretion disk models in which the assumption of Keplerian rotation law is relaxed and the 
radial advection term in the energy equation is included explicitly, it has been revealed 
that advection becomes a dominant cooling process when the mass accretion rate is larger 
than the Eddington limit (see, e.g., Abramowicz et al 1988; Kato, Honma, & Matsumoto 
1988; Chen & Taam 1993). Under that situation, the disk is both thermally and viscously 
stable even though the disk is radiation pressure dominated. 

In optically thin hot accretion disk models, it has long been assumed that advection is 
not important and was neglected. Optically thin hot accretion disks have been proposed to 
explain the hard spectrum observed in black hole candidates such as Cyg X-l. Early work 
was published in mid-1970s (see Shapiro, Lightman, & Eardley 1976; Liang & Price 1977; 
and review by Eardley et al 1978), and these concentrate on the spectral aspects. Recent 
work includes Kusunose & Takahara (1988, 1989, 1990), Tritz & Tsuruta (1989), White & 
Lightman (1989, 1990), Bjornsson & Svensson (1991a, b, 1992), and Kusunose & Mineshige 
(1992). These studies are focused on the microphysical processes such as electron-positron 
pair production and annihilation of the hot plasma (see, e.g., Svensson 1982, 1984). In all 
these models, global heat transport due to radial advection has not been included. Disks 
with advection have been studied recently by Narayan & Popham (1993), Narayan & Yi 
(1994) and Abramowicz et al (1995). Narayan and Popham (1993) considered accretion 
on to white dwarfs and got advection dominated solutions for low accretion rates when the 
local cooling is inefficient. Abramowicz et al (1995) obtained thermal equilibrium curves of 
accretion disks at a given radius with a prior assumed advection cooling term. They showed 
that an advection dominated, thermally stable, solution co-exists with the local cooling 
dominated thermally unstable hot disk solution. These results may be important not only 
on the theoretical aspect but also for the understanding of the observed phenomena from 
accretion powered sources. 







Accretion disks with advection included are governed by a set of stiff differential 
equations. They are numerically very difficult to solve as is proven in the case of optically 
thick slim accretion disks surrounding black holes (see, e.g., Abramowicz et al 1988; Kato, 
Honma, & Matsumoto 1988; Chen & Taam 1993) and boundary layers (e.g., Narayan & 
Popham 1993). In the approach of Abramowicz et al (1995), a Keplerian disk is assumed 
and the advection term is localized by assuming the dimensionless radial differential factor, 
£ (see eq. [10]), to be unity. Thus local (i.e., fixed radius) solutions of the surface density 
can be calculated for given mass accretion rates with one algebraic equation. However, 
inconsistency may occur since the advection is not calculated from the global structure 
of the disk which has not been constructed. In this paper, we solve the global energy 
conservation equation explicitly and so the consistency in regard to the radial advection 
heat transport is insured. In the next section, we briefly describe the assumptions and 
formulation used in the study. In §3, we investigate the role of advection quantatively. 
For advection dominated solutions, an analytical formula will be derived in a self-similar 
approximation. Detailed numerical results on the radial structure of the accretion disks 
will be presented. Discussion is presented in the last section. 



2. FORMULATION 

We consider an accretion disk which is axisymmetric, non self-gravitating, opti- 
cally thin and geometrically slim so that it can be described by the vertically inte- 
grated equations. In this approximation, the Keplerian rotational angular velocity is 



fife = \J GM I R(R — Rq) 2 under the pseudo-Newtonian potential $ = — GM/(R — Rg), 
where M is the mass of the central object and Rq = 2GM/c 2 is the Schwarzschild ra- 
dius (Paczyhski & Wiita 1980). We assume that the angular velocity can be expressed as 
O = \J GM I R 3 /u, thus for Keplerian disks, u = 1 — Rq/R. In a steady state, the angular 
momentum equation is reduced to 

fV = ^fg-\ (i) 

where M, E and v are the mass accretion rate, the surface density, and the kinematic viscos- 
ity respectively. The correction terms are / = 1 — 9fi*/(fir 2 ) and g = — (2/3)(<fln fi/c/ln R). 
Here fi* = Q,(3Rq) and r = R/Rq. We use a standard a-model viscosity prescription 
(Shakura & Sunyaev 1973): 

H, (2) 
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v = —ac 
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where a is a constant, c s = \fp~fp is the local sound speed, and H = c s /fifc is the half- 
thickness of the disk. Here p and p are the density and total pressure of the disk re- 
spectively, E = 2Hp. In a steady state, the mass conservation equation can be written 

as 

M = -2vri?Eu r , (3) 

where v r is the radial velocity. The steady state energy conservation equation is repre- 
sented by the balance between the local viscous heating, Q+, the local radiative cooling 
in the vertical direction, Q- } and the global heat transport (radial advection), Q a dv It is 
expressed as 

Q+ = Q- + Qadv (4) 
The viscous heating rate per unit area is given by the standard formula 

3G MM _,_ 

The advection cooling rate is taken in a form (see, e.g. Chen & Taam 1993): 



n - v £ 

P 



4-3/3dlnT m dlnS 

r, -l dR ~ {A ~ m ^R~ 



(6) 



where T is the mid-plane temperature of the disk, (3 = p g /p is the ratio of the gas pressure 
to the total pressure, T3 = 1 + (4 — 3/3)(7 — l)/[(3 + 12(7 ~~ -Q(/^ ~~ an d 7 1S the ratio of 
specific heats. For optically thin disks, p = p g = IZpT '/ p and (3 = 1, where p is the mean 
molecular weight assumed to be 0.617. We shall assume that the local radiative cooling is 
provided by optically thin thermal bremmstrahlung with emissivity (erg s _1 cm -2 ), 

Q- = Qbre m = 1-24 x 10 21 Hp 2 T 1 / 2 A } (7) 

where A > 1 is the Compton enhancement factor. Combining the above equations and the 
vertical hydrostatic equilibrium relations we obtain (see also Abramowicz et al 1995), 

T = 1.39 x 10 13 (aS)- 1 mr- 3/2 /fl'" 1 ?" 1 , (8) 
im 2 - 0.361r 1/2 ^ 2 gcu" 2 (aE)m + 3.462 x 10" 6 Aa" 2 r 2 /" 1 ^(aE) 3 = 0, (9) 

and 

1 dlnT dlnE 
^ ~ ~ 7 - 1 dlni? + dlniT 







Here q = 1 — Rq / R } rh = M/Me } and Me = AttGM / (cn es ) is defined as the Eddington 
rate, where n es = 0.34. 

3. STRUCTURE OF THE DISK 

For a given mass accretion rate m, the radial disk structure (i.e., the surface density 
and temperature) can be obtained numerically from equations (8)-(10). However, the 
characteristic properties of the solution can be revealed under some approximations. In 
equation (9), one notes that the first, second and third terms correspond to the radial 
advection heat transport, the viscous heating and the local radiative cooling respectively. 
Thus by neglecting the first term, the usual purely local radiative cooling hot accretion 
disk solution can be recovered. Furthermore, equation (9) alone can be solved by assuming 
a constant £ (Abramowicz et al 1995). For a given mass accretion rate rh and a fixed radius 
r, this equation is a cubic equation in E and in general has three roots. They are either 
three real roots or one real and two complex conjugates. In the case of a postive £, which 
corresponds to a cooling radial advection (i.e., to a negative entropy gradient), a critical 
mass accretion rate exists (see Abramowicz et al 1995): 

m max = 1.7 x 10 3 a 2 r- 1 / 2 A- 1 fg 5 q 3 u- 6 C\ (11) 

which decreases for a large Compton enhancement factor. For rh < m max there are one 
unphysical negative root of E and two postive roots corresponding to the local cooling and 
advection cooling dominated solutions. For rh > m max , however, there are one negative 
root and two complex conjugates. They are all unphysical. In the case of a negative £, 
which corresponds to a heating radial advection, there is only one physical positive solution 
and the other two are negative and thus unphysical. 

It should be noted that, the mass of the central object does not appear in equations 
(8)-(10). Thus, unlike optically thick accretion disks, the temperature and surface density 
of hot optically thin accretion disk do not depend on the mass of the central object if 
they are expressed in terms of R/Rq and M/Me- This occurs because the disk is gas 
pressure dominated, thus (H J R) 2 oc T and no M is explicitly involved. Since the heating 
rate, Q+, scales to M _1 , it requires the cooling rates (both the local and the advection) 
be proportional to M _1 also. This is the case for advection, Q a dv oc M _1 (H J R) 2 , and 
for bremmstrahlung, Qbrem on M _1 (H/R) -1 , but not for optically thick cooling in the 
radiative diffusion approximation. 
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3.1. The Role of Advection 

We concentrate on the case where the radial advection is a cooling process and the 
mass accretion rate is less than the critical value so that two solutions can exist. To see 
the effect of the radial advection in each of these two solutions, the Q a dv/Q+ ~ M relation 
is shown in Figure la. For simplicity, we have assumed that A = 1, u = 1 — Rq / R } £ = 1, 
a = 0.1, and r = 10. It is seen from the right vertical line (which corresponds to the upper- 
left branch line of the M(E) relation in Fig. lb) that the disk is fully advection dominated 
for a wide range of mass accretion rates, that is, Q a dv/Q+ ^ 0.999. For mass accretion 
rates close to m max , the advective and local cooling processes are equally important (see 
the horizotal line which corresponds to the turning point and a small section of the lower 
branch of the M(E) relation). However, for the rest of the solutions on the lower branch of 
M(E) (Fig. lb), the advection is negligible, that is, Q a dv/Q+ ~ 0.001 (see the left vertical 
line). This result suggests that, for mass accretion rates less than m, max , the usual hot 
accretion disk solutions can be constructed in a very good approximation by neglecting the 
radial advection heat transport effect, although they are thermally unstable. On the other 
hand, for the advection dominated solutions, the detailed local cooling process will not 
affect the disk structure very much. It can be easily checked that above the advection line, 
Qadv > Q+i whereas below this line, Q a dv < Q+- Thus, for a perturbation of M (increase 
M but keep E unchanged), or effectively an increase of temperature, the cooling rate will 
exceed the heating rate and so the gas cools down and the original thermal equilibrium is 
restored. In other words, the advection dominated solution is thermally stable. 

3.2. Self-similar Advection Solution 

The dimensionless radial differential factor, £, can be estimated from a self-similar 
solution of accreion disks, which can be obtained by assuming the correction terms, /, g, 
q } and u to be constants independent of the radius of the disk. Self-similar solutions of 
accretion disks have been studied by various researchers (see, e.g., Begelman & Meier 1982; 
Spruit et al 1987; Liang 1988; Narayan & Yi 1994). Although self-similar solution has its 
own limitations, such as that it cannot be applied to the inner regions of the disk due to the 
restrictions of the boundary, it has the advantange of reducing the global radial advection 
to a localized quantity. The self-similar solution for advection dominated accretion disks 
can be obtained in the following procedure (see also Narayan & Yi 1994). 

Under the assumption of the correction terms being constants, from equation (8) we 



have, TE oc r _3//2 . From equation (9) by neglecting the third local cooling term, we have 
E oc r -1//2 . Thus, T oc r _1 and 

C = — -, (12) 

^ 7-1 2' v ; 

which is independent on the disk parameters such as the mass accretion rate, the mass of 
the central object, and the viscous parameter a. It gives £ = 1 for 7 = 5/3. It should be 
noted that the advection term we adopted is vertically integrated, and uncertainties are 
unavoidably introduced (see Abramowicz et al 1995 and Narayan & Yi 1994 for a different 
definition of £). However, this result can be regarded as a clear demonstration that radial 
advection in accretion disks is a cooling process under the self-similar assumption. 



3.3. The Radial Disk Structure 

We have solved equations (8)-(9) with an iteration technique for the local radiative 
cooling dominated solutions. The calculation starts at the out edge of the disk and con- 
tinues inwards. By first assuming a numerical value for £, we calculate two set of solutions 
of E and T at radii r + Ar/2 and r — Ar/2 from equations (8) and (9). A new £ is then 
obtained from equation (10). With this new value of £, we repeat the calculation of E and 
T. Finally a converged £ will be found. At the next radius, the previous £ is used as the 
first guess. 

For the advection dominated solutions, we use a shooting method to integrate the 
ordinary differential equation of E derived from equations (8)-(9). The initial starting point 
for the integration is located at a very large distance from the central object. However, 
solutions in regions, for example, inside 1000i?G, which we considered here, is not sensitive 
to the exact location of the starting point. The initial starting value of E is obtained 
through iterations as stated above. For definiteness, we assume here that the disk is a 
Keplerian one. 

We first present the relation of the dimensionless radial differential factor, £, with the 
radius of the disk and the mass accretion rate (Fig. 2). We see that for both type of 
solutions, namely, the advection dominated (solid lines) and the local cooling dominated 
(dotted lines), £ does not depend on the mass accretion rate very much; there is only 
a negligible decrease for a wide range of rh (from 0.1 to 0.001). This result confirms 
the assumption made by Abramowicz et al (1995) that £ is a constant on the thermal 
equilibrium curve of M(E) for a given radius of the disk (see Figure lb, the lower-upper 
branch shaped curve). It implies that the shape of the equilibrium curve will not change in 
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more detailed models. One also notices that £ of advection dominated solution are larger 
than that of local cooling dominated solution and it approaches to the self-similar solution, 
£ = 1, at larger radii of the disk. For the inner regions of the disk, £ decreases, however, 
due to the zero torque condition applied at the inner boundary of the disk. 

The temperature of advection dominated accretion disks is much hotter than that of 
local cooling dominated ones and it does not depend on the mass accretion rate (Fig. 3). 
This reflects the fact that in this case the disk heats up to the virial temperature (Narayan 
& Popham 1993). On the other hand, the surface density of advection dominated accretion 
disks is much smaller than that of local cooling dominated disks and it depends on the mass 
accretion rate more sensitively(Fig. 4). For advection dominated solutions, the surface 
density will not depend on the local radiative process much since it is determined primarily 
by the radial advection and the viscous heating terms and the local cooling term has a 
negligible effect (see eq. [9]). The weak dependence of the temperature on the mass 
accretion rate and the viscous parameter, a, is obvious from equations (8) and (9) for 
advection dominated disks, in which 

aS = 2.773mr- 1/2 g- 2 q- 1 co 2 C, (13) 

and 

T = 5.02 x lO^r" 1 /^" 2 ^" 1 - (14) 

It also shows that, taking a constant £ = 1 in the whole disk will give lower temperatures 
in the innermost regions as one then has less energy to heat the gas due to an overestimate 
of the advection. In the local radiative cooling case, on the other hand, we have 

S = 3.23 x lOW/V-^f/y/y 72 ^- 1 ^ 1 / 2 , (15) 

and 

aT = 4.30 x lO 10 ^ 1 /^- 3 / 4 / 172 ^" 372 ?" 372 ^ 172 , (16) 

which shows instead that the temperature depends on both the viscosity parameter and 
the mass accretion rate. 

The extremely high temperature of the advection dominated accretion disks infers the 
presence of various microphysical processes in the hot plasma. The equations governing 
theses processes depend on the assumptions of the model and may include such as the 
electron-positron pair balance, the photon blance, the pair energy balance, and the proton 
energy balance. In the case of a spherical homogeneous hot plasma cloud (HPC), two 
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parameters can be used conveniently and they fully determine the physical conditions 
in the HPC (see, e.g., Svensson 1982, 1984; Bjornsson & Svensson 1991b). These two 
parameters are defined as the proton optical depth, t p , and the global compactness, /: 

t p = cr T n p H; I = (17) 

where gt = 6.65 X 10 _25 cm 2 is the Thomoson cross section, n p = p/m p is the proton 
number density, m e is the electron mass, m p is the proton mass, and H and L are the size 
and the luminosity of the HPC respectively. In the case of an accretion disk, Bjornsson 
& Svensson (1991b) identified H with the local scale height of the disk and replaced the 
global HPC compactness parameter with a local disk compactness parameter defined as 



l=^^- 3 H*Q-. (18) 



From equations (7) and (8), we then have, 

t p = 0.2S, (19) 

and 

1 = 3.948 x 10 4 a _1 SmmA/s _1 . (20) 
With equations (8) and (19), we can express the Comptonization parameter in E and m, 

y c = 40 Max(rp,r 2 ) = 1.875 x 10 3 a _1 m Max(l, 0.2S) r _3/2 fg' 1 q' 1 , (21) 

where = kT / (m e c 2 ) is the dimensionless temperature. 

The loci of / = 10, t p = 10 -2 , and y c = 1 along with the thermal equilibra are shown 
on the E ~ M plane in Figure 5 for three different viscosity parameters, a = 0.001, 0.01 
and 0.1 at radius r = 10. Since we have expressed M in unit of «Me, these loci and 
also the advection dominated branch do not change with a (see eqs. [13], [19-21]). If 
we assume electron-positron pair production is important for the compactness parameter, 
/ > 10 (e.g., Svensson 1984), then this condition is M > 0.01«Me for advection dominated 
solutions, and which is higher than that of the local cooling dominated solutions. On the 
other hand, for both kinds of solutions, Comptonization effect becomes important only for 
M > O.IckMe. The locus of y c = 1 becomes flat (i.e., independent of E) when the proton 
optical depth, t p , is smaller than unity. 

To examine the importance of these processes at different regions of the disk, the 
variations of /, t p , y c , and 0, with respect to the radius of the disk are shown in Figure 6 
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for three different mass accretion rates, m = 0.001, 0.01, and 0.1, with viscosity parameter, 
a = 0.1. Here the disk is advection dominated. It is seen that the disk compactness 
parameter does not change very much with the radius of the disk and so the electron- 
positron pair process is important in the entire advection dominated accretion disk for 
M > O.OIckMe. The Comptonization parameter, y c , however, decreases rapidly with 
increasing disk radius. This is due to the fact that both the temperature, 0, and the 
proton optical depth, T p} decrease with radius. This result may imply that the importance 
of Comptonization is restricted to the inner regions of the disk. Finally, we also see that 
relativistic effects may be important since > 1 (e.g. Svensson 1984). 

Advection dominated accretion disks are most probably in the two-temperature regime 
since > 1 (e.g. Svensson 1984). For a two-temperature accretion disk, the temperature 
shown here corresponds approximately to the sum of the proton (ion) temperature, T p} and 
the electron temperature, T e . Since T p > T e , we have T T p . The electron temperature 
can be determined, for example, through the proton-electron Coulomb coupling rate (see 
e.g., Guilbert & Stepney 1985), 

F = 3H (IT te) aTC{ln mT ' - T - ] te) 2 {2z + (22) 

where z = n+/n p} n_|_ is the positron number density and In A 20 is the Coulomb 
logarithm. A simple estimate can show that the electron temperature is approximately 
one order of magnitude smaller than the proton temperature which is in the range of 10 n K 
to 10 12 K. 

4. DISCUSSION 

The global radial structure of hot optically thin accretion disks with advection in- 
cluded explicitly has been calculated. It is shown that advection is generally a cooling 
process which corresponds to a negative radial entropy gradient. For a given mass accre- 
tion rate less than a critical value, there are two solutions in which the advection can be 
either a dominant cooling process or negligible. Whereas for mass accretion rates above 
the critical value no hot optically thin accretion disk solution can exist. This result does 
not depend on the detailed local radiative cooling process. However, it requires that the 
purely local cooling solution has an approximately straight line relation of M(E) (at a 
fixed radius on the logarithmic plane) with a slope steeper than d( log M)/d( log E) = 1/2 
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or c/(log M)/d( log E) < 0. In fact this is the case for non-pair bremsstrahhmg and Comp- 
tonized soft photon cooling solutions (see, e.g., Wandel & Liang 1991, Luo & Liang 1994), 
and for electron-positron pair accretion disk solutions (see, e.g. Bjornsson & Svensson 
1992; Kusunose & Mineshige 1992). 

We have used bremsstrahlung radiation mechanism for the local cooling process. At 
temperature exceeds 10 9 K, other radiation mechanisms can often cool the electrons much 
more efficiently (see e.g., Rees et al 1982). This will result the local cooling dominated 
lower-branch of the thermal equilibria to move in the up-left direction and the critical mass 
accretion rate, m max , will decrease (see Fig. 5). However, these local cooling processes 
will not affect the advection dominated upper-branch much. This is because, as have been 
demonstrated in §3.1, that for advection dominated solutions, the local cooling process 
has a negligible effect (Q a dv/Q+ ~ 0.999). The detailed study of the microphysics of the 
hot plasma is beyond the scope of the present paper. However, we have suggested that 
the disk is in the two-temperature regime. The temperature calculated in this paper is 
more closely related to the proton instead of the electron temperature. A self-consistent 
electron temperature can be obtained by considering the detailed radiative microphysics 
of the hot plasma. We have also shown that complicated radiation processes are present 
most probably in accretion disks with large viscosity parameters and high mass accretion 
rates (see Fig. 5). Thus, our results are more applicable to accretion disks with small 
viscosity parameters and low mass accretion rates. 

The advection dominated accretion disks can hardly be regarded as in a geometrically 
thin configuration. By considering two-dimensional (vertical and radial) effects, results 
obtained here will probably be changed in some way. Furthermore, the non-Keplerian effect 
which is obviously important due to the high local sound speed, and the transonic nature 
of accretion flows near the vicinity of the black hole will also play roles in determining the 
global disk structure. All these will be the subjects of our future investigations. 

The author thanks Profs. Marek Abramowicz, Jean-Pierre Lasota and Ronald Taam for 
stimulating discussions and useful suggestions. Comments from them and Drs. Gunnlau- 
gur Bjornsson and Fredrik Wallinder are highly appreciated. 
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FIGURE CAPTIONS 

Figure 1. — a) The variation of Q a dv/Q+ with respect to rh. Here, A = 1, u = 1 — Rq / -R, 
£ = 1, a = 0.1, and r = 10. b) The m(E) relation. Note that on the upper branch, the 
advection is fully dominated whereas on the lower branch, the advection is negligible. 

Figure 2. — The variation of the dimensionless advection factor £ with respect to the 
radius of the disk and the mass accretion rate. The solid and dotted lines are for advection 
and local cooling dominated respectively. Here, a = 0.1 and rh = 0.001 — 0.1. Note that £ 
does not depend on the mass accretion rate very much. There is only a negligible decrease 
for a wide range of rh (0.1 to 0.001). 

Figure 3. — The radial temperature structure of hot optically thin accretion disks. The 
dotted lines from top down represent the local cooling dominated solutions for rh = 0.1, 0.01 
and 0.001 respectively. The solid lines represent the advection dominated solutions which 
do not depend on the mass accretion rate. 

Figure 4. — The radial surface density distribution of hot optically thin accretion disks. 
The dotted lines from top down represent the local cooling dominated solutions for rh = 
0.1,0.01 and 0.001 respectively. The solid lines from top down represent the advection 
dominated solutions for rh = 0.1,0.01 and 0.001 respectively. 

Figure 5. — The loci of / = 10, t p = 10 -2 , and y c = 1 along with the thermal equilibra. 
Here, r = 10 is assumed. Note that pair production is important for M > O.OIckMe 
(advection dominated solutions); and Comptonization effect becomes important for M > 
O.IckMe (both kinds of solutions). 

Figure 6. — The variations of /, t p , y c , and 0, with respect to the radius of the disk. Here 
the disk is advection dominated and a = 0.1. The solid, dotted, and dashed lines represent 
rh = 0.1, 0.01, and 0.001 respectively. Note that the disk compactness parameter does not 
change very much with the radius while the Comptonization parameter decreases rapidly. 



